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Abstract 

Removal of the quenched approximation in the mechanism which produced an analytic estimate 
of quark-binding potentials, along with a reasonable conjecture of the color structure of the nucleon 



r^ . formed by such a binding potential, is shown to generate an effective, nucleon scattering and 

Ph. 
f->i ' binding potential. The mass-scale factor on the order of the pion mass, previously introduced 

D ■ 
,^ ' to define transverse imprecision of quark coordinates, is again used, while the strength of the 



potential is proportional to the square of a renormalized QCD coupling constant. The potential 
so derived does not include corrections due to spin, angular momentum, nucleon structure, and 
electroweak interactions; rather, it is qualitative in nature, showing how Nuclear Physics can arise 
from fundamental QCD. 
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I. INTRODUCTION 

In previous papers [lH5|] a new, analytic, non-perturbative, gauge-invariant approach to 
QCD has been defined and used to give a simple estimation of a quark-binding poten- 
tial, V{r) ~ ^/i (r /i)(^+^\ where ^ is a small, real, positive parameter, ^ <^ 1, introduced 
phenomenologically into that function which guarantees transverse quark-coordinate impre- 
cision; ^ and fi are to be determined in terms of the pion and nucleon masses, as noted in 
Ref. [5|; /i is a mass scale parameter understood to be on the order of the pion mass, m.,,. 

In this paper, we begin with the concept of three bound quarks scattering against another 
triad of three bound quarks, with their full, non-perturbative exchanges of gluons between 
all the quarks taking place. We assume that those triads that are initially bound remain 
bound at all times, which carries the implication that the multiple gluons exchanged between 
these nucleons do not change the overall color-singlet nature of each nucleon. We neglect 
all electroweak interactions; and to further simplify the analysis assume that this nucleon 
scattering takes place at high relative velocities, so that a simplifying eikonal description of 
the scattering may be used. We further simplify the analysis by neglecting spin effects - 
which could be inserted if desired - and aim for a simple, qualitative picture of how forces 
between nucleons can arise, starting from the basic fundamentals of QCD. 

In one sense, however, our eikonal model must be made more complicated than those 
quenched models used previously, for it turns out that one must here retain at least the 
simplest effects of the closed-quark-loop, or vacuum functional, L[yl]. The Physics underlying 
this requirement follows because the forces which arise between quarks due to the multiple 
exchange of gluons are strong, tending to bind, for impact parameter b separations on 
the order of l//i, but they fall off rapidly with increasing b > 1/yU. (At distances b <^ 
1/yU large color fluctuations tend to reduce the value of any amplitude, and that could 
very well correspond to a non-perturbative translation of asymptotic freedom). How can 
nucleons, whose internal structures are deflned at distances b ~ 1/yU, then feel strong forces 
at separations b > l//i? 

The answer is that vacuum loops, deflned by L[yl], can stretch in the transverse directions, 
and can serve to transmit the multiple gluon interactions across larger values of impact 
parameter; a gluon "bundle" from one nucleon attaches itself to one point on the loop, while 
another bundle of gluons passes from a second point on the loop (at a signiflcant transverse 



distance from the first) to the other nucleon. Although this passage of momentum via a 
closed vacuum loop changes the interaction somewhat, an essential "short-range" interaction 
is produced at distances larger than would be possible by quark-to-quark passage alone. 
And if, as this loop is stretched in a transverse direction, one transverse side of the loop 
corresponds to a quark and the other to an antiquark, one has the image of an effective pion 
being exchanged between scattering nucleons. 

Effective Locality (EL) is the reason that multiple gluon exchange appears to leave or 
arrive at a single space-time point on a "quark line", in the "bundle- diagrams" to follow, 
which correspond to the sums over infinite numbers of conventional Feynman graphs. Such 
sums are possible because they can be represented by the equivalent of Gaussian functional 
integration over products of factors of Gc[v4] and L[y4], where Gc[v4] represents a quark 
propagating in a "classical" gluon vector potential y4"(x); and that integration can be car- 
ried out exactly because there exist Fradkin representations of these two functionals which 
are Gaussian in their A-dependence. Results are then expressed in terms of these Frad- 
kin representations, which are just Potential Theory constructs, and have relatively simple 
approximations in different physical situations. In this paper, for reasons of simplicity, we 
shall replace the Gc[A] of each quark in a nucleon, and the single Gd-A] which models that 
nucleon, by its high-energy eikonal limit, and then connect the gluon bundles emitted by 
each nucleon to two, and only two points on a single loop. More complicated loop structures 
are certainly possible, and should be investigated, but this is the simplest representation of 
"effective pion exchange between nucleons". 

Finally, mention should be made of the relative simplicity of this approach, compared to 
other well-known and long-studied methods of calculation in QCD. Again, it is the unex- 
pected EL, appearing automatically after the non-perturbative sum over all possible, gauge- 
invariant, gluon exchanges between quarks has occurred, which is responsible for the huge 
simplifications obtained as Halpern's functional integrals are here reduced to a finite set of 
ordinary integrals |6| amenable to computer evaluation, and here estimated in the simplest 
way possible. One has long believed in the Principle of " Conservation of Difficulty''' , when 
calculating higher-order effects in QED or any Abelian theory; but for non-Abelian field 
theories, approached in the manner we suggest, that Principle is not true. 



II. FORMULATION 

We begin by first presenting a bird's eye view of tfie detailed calculations that follow 
immediately, expressed in terms of the analysis of |2|, and in particular to the discussion 
centering about its Eqs. (20) - (23), where the color coordinates of each of the two Gc[A], 
each representing a scattering quark or antiquark, are discussed, in connection with the 
eikonal forms at high energies, taken by the exact Fradkin representation of that Green's 
function. Now generalize that process to the scattering of a triad of quarks bound into a 
singlet nucleon with another such nucleon. For each nucleon, there will occur the product of 
three such terms as written in Eq. (23) of |2| , each with the same fraction of that nucleon's 
momentum, and with their color weightings Q°' restricted so as to insure that the three 
bound quarks comprising each nucleon remain in a color-singlet state. 

We replace the description of that combination by that of a nucleon, of momentum p 
and effective color weighting Q°-, so defined such that only those combinations of Gell-Mann 
matrices of each of the basic quarks corresponds to gluons which may be exchanged and 
so preserve each triad of quarks in its bound, color-singlet state. Note that all such color- 
singlet exchanges can be absorbed and emitted by the quark line comprising the loop L[y4], 
for the Fradkin representation of L[A] contains a trace over all possible combinations of 
color coordinates. By this simplification, we replace the essence of a 6-body problem by 
a 2-body problem; and this has the consequence that our subsequent estimation of the 
nucleon-nucleon binding potential produces a qualitative description of how nuclear forces 
can arise from basic QCD. 

We next refer the reader to Section IV of Ref. J3j, and in particular to the func- 
tional operations of its Eq. (30), with attention drawn to the translation operator 
exp[— J Q ■ {gf ■ x)^^ ■ ^] acting upon exp {L[A]}. Here, Q refers to the coordinates of 
the two nucleons, and the translation operator inserts that dependence in a well-defined 
way into the i^[A] written in Eq. (All) of that paper. The remaining linkage operator of 



(30) generates a functional cluster expansion, discussed and derived in Ref. [13| and |l4| : 
for our purposes, involving but a single loop and suppressing any gluon bundles exchanged 
across that loop, this remaining linkage operation can be neglected. This simplified analysis, 
which will require a simple renormalization, is sufficient to produce a reasonable, qualitative 
nuclear potential from basic 'realistic' QCD. By 'realistic', it meant a formulation of QCD 



which contains from its inception that asymptotic quarks and/or antiquarks are found only 
in bound states, and hence their transverse coordinates cannot be specified exactly. In this 
paper, such transverse imprecision follows from the defining arguments found in [4|. 

The specific steps of this analysis follow from a statement of the Generating Functional 
(GF) derived in ^, 
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which is nothing else but a particular, and gauge-invariant rearrangement of the 
Schwinger/Symanzik GF derived more than a half-century ago, and applied to QCD. (Such 
a representation of gauge invariance is dependent upon cubic and quartic gluon interac- 
tions, and is not possible for QED.) Here, the covariant, causal gluon propagator is given, 
in momentum space, by 
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with C an arbitrary gauge parameter, while exp {D^i} represents the "linkage operator", 
with 
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This functional differentiation formalism has the great advantage over the often more cus- 
tomary Functional Integral (FI) over gluon fluctuations, in that there is no need for concern 



about spurious (Gribov) replication of gauge copies; it is gauge- invariant |l8 



from the very 



beginning, and made so by means of a small observation overlooked for decades, as described 
in [21. 

The functional Gc[v4] = [m + 7^ (9^ — z^fA^jA")] represents a quark propagating in 
the presence of a specified, "classical" field Aj^(a;), while L[A] = Tr In [1 — igi^j ■ A ■ X) Sc], 
Sc = Gc[0], denotes the "vacuum" functional; both functionals have Fradkin representations 
essentially Gaussian in their A- dependence, which permits the linkage operations needed for 
the specific amplitudes desired (obtained by functional differentiation with respect to the 
sources ri{x), f]{y), j'^{z)) to be obtained exactly; and this, in turn, corresponds to the 
summation of an infinite number of classes of Feynman graphs, each class containing an 
infinite number of graphs, with the results expressed in terms of the parameters of the 



Fradkin representations for Gc and L. In this paper, for simplicity, we shall replace the 
Fradkin representation for Gc[v4] by its high-energy, eikonal approximation; and, again for 
simplicity, restrict consideration of L[A] to its simplest form in which quark spins have been 
suppressed. Finally, as discussed and derived in [4|, the "transverse averaging" needed for 
"realistic" QCD, after the non-perturbative sums over Feynman graphs have been computed, 
can be very simply inserted as an intermediate step by the replacement 

A;(x - u{s')) -^ / d^sl a{x± - xl) A^^ix' - u{s')), x'^ = (zxq; fl, xi^), (4) 

where Ci{x±—x'^ is a real, symmetric function, expressing the impossibility of defining precise 
transverse coordinates of bound quarks and/or antiquarks. (Originally, as in |^, only the 
quark color current operator was endowed with such transverse imprecision; but because 
of the assumed symmetry of the a(x_L — a;'j_), in that part of the Action operator coupling 
such color current to gluons, it is permissible to interchange the roles of transversally-shifted 
fields, and imagine that it is the coordinate dependence in A which has been shifted. In 
reality, no such shift of the F" (x) have been performed.) 

It may also be noted that cubic and quartic gluon interactions are fully included in this 
formalism, and are represented by the Halpern FI over d[x]- And because of the remark- 
able property of EL, alluded to above, which appears after the non-perturbative sums are 
performed, the Halpern FI of ([H) reduces to a finite set of ordinary integrals 6|, which are 
amenable to numerical computation. In this paper, again for simplicity, we estimate these 
integrals by means of an approximate Gaussian evaluation. 

We now ask the reader to imagine that functional derivatives are taken with respect to six 
pairs of r^, ff sources, "bringing down" six G^i'^^Kv \A)^ which we divide into two groups 
of three, 

G,(x«,y;^V) Gc{xf\yf\A) G,{xf\yf\A) (5) 

for nucleon I, and a similar triad with coordinate superscripts (4), (5), and (6) for nucleon 
E, beginning the computation as if we were calculating a six-quark scattering amplitude. 
Each Gc[v4] will bring to its triad the A-dependence associated with an eikonal/high-energy 
limit of its exact Fradkin representation, of form 



exp 






(6) 



We now introduce the bound-state nature of each triad of quarks by first suppressing the 
coordinate superscripts, ?/^J^ — )■ 2/(1 j), and Pn\\ — ?• (1/3) P/u), since each quark must have 
essentially the same space-time and momentum coordinates if its nucleon is to remain intact. 
This means that the product of the three factors of Eq. ([H]) which are now a property of 
each nucleon may be written as 



exp 
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(7) 



We emphasize that ([7]) refers to the A-dependence of each nucleon after the linkage 
operations binding each triad of quarks have been performed, as discussed in [5|, which 
analysis we here suppress. Further, for each nucleon to remain bound for all (proper) times, 
there must exist a relation between the f2"''^*''(s) such that only color singlets are exchanged 
between nucleons I and E, and this relation should be independent of s. Since the fl'^'^'^^ 
define the Gell-Mann matrices A"'*-*-* inserted between initial and final nucleon states, there 
must be a relation between the fi"'*^*-' guaranteeing that each nucleon remains a color singlet. 
We thus simplify Eq. ([7]) by introducing Q?j J^^ 
Eq. (I7l) becomes 

Pm,(i,i) 
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X]j=i ^aV)' ^^^ re-scaling s — )■ 3s, so that 
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A modification, representing the "realistic" QCD defined and used in the two preceding 
|4, |5|, replaces in the exact Fradkin representation each v4"(?/' — u(s')) by 



papers 



d'y'^a{y^-y'^)Al{y'-u{s')), 



(9) 



where y'^ = (i^/g; yj_, y^) and y^ = {iy^; y±, y^) represents the coordinate of an individual nu- 
cleon. Since we are assuming strict binding of each nucleon, the only transverse imprecision 
we must specify for this analysis is that between the quarks of one nucleon and those of the 
other; and since we have assumed that such gluon exchanges are not in any way to disrupt 
the binding of quarks within each nucleon, we shall invoke transverse imprecision for the y'^ 
coordinates of A'^{y' — u{s')) of (jH]), replacing (^ by 



exp 



d'«^2;:,(i,n)(^)A^(«;) 



(10) 



where 



2^,(1,1) H= M^)^(W (11) 
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In so doing, we have replaced the Fradkin coordinates u^^^s) by the eikonal combinations — 
or ^^, and have neglected quark spin dependence. 

The binding process effectively transforms each triad of quark Green's functions into 
a single nucleon Green's function of mass M and 4-momentum p^ where the exponential 
of (fTOj) . the effective 'rehc' of its original three quarks, is retained for subsequent use in 
calculating the interaction between both nucleons. Each nucleon's Green's function now 
contributes to the eikonal scattering amplitude of the two nucleons, and in their CM takes 



on the standard form 19 1 



r(.,t) = ^yd^e-^[l-e^^(M], (12) 

where g is the momentum transfer of this scattering process, s = —{pj + pj)^, t = — (pj — 
p'^y = — g^, and where 

^)= A/" fd[x]e-i^^' (13) 



^^I 



A->-0 



with normalization constant A/" defined such that X — )■ for (7 — )■ 0. 
The linkage operation of (IT31) then has the Gaussian form 



e-U^-^^-TA .gt/A■/C.A+^/7^.AgL[A] 
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(14) 



where JC"^^ = gf'^^x'L + (D^ ^)ll and W^ = d^x% - Qi,^ - Qi,^- As in the passage from 
(21) to (31) of Ref. (3|, the functional operation may be carried through exactly, yielding 
for dH]), 



A^O 



The first line of (llSp may be rewritten as 

[det {gf ■ xT^ ■ exp L'- j [dx -Qi- Qi] ■ {gf ■ xT' ■ \dx - Qi - Qi]}, (16) 



and we here rely on the strong couphng hmit of (? ^ 1 , keeping only the terms Qi^i (this is 
really not necessary, but it simplifies the analysis; if the (9x-terms are retained, the normal- 
ization integrals become more complicated, but the thrust of the procedure is the same). 

Furthermore, the terms of fITB]) proportional to two factors of Qi and to two factors of 
Qn are "self-energy" corrections to the respective nucleon propagators, and they will be 
suppressed, since we are here interested only in the interaction of one nucleon upon the 
other. A similar remark may be made for those terms containing a single factor of dx and 
either Qi or Qn, for they correspond to " tadpole" -hke structures attached to either nucleon, 
and are not relevant here. With these simplifications, (ITBl) is replaced by 



[det {gf -x)] ^ ■ exp 



Qi-igf-xr-QiL 



(17) 



which, except for different color factors, has the form of the eikonal function describing the 
interaction between a pair of quarks. 

For the impact parameter range between nucleons in which we are interested, it turns out 
that fITTj) gives an unimportant contribution to the nucleon-nucleon potential; for simplicity, 
we here neglect it, in contrast to the true source of that potential, which arises from the 
action of the linkage/displacement operators of flTSl) upon exp {L[y4]}. It is worth mentioning 
that this makes for an important difference with the quark binding potential evaluated in 
Ref. y. 

Denoting the linkage operator of ( 1T5|) by 



exp [^a] = exp 



"2 / ^^•( ^/-X) ■ ^^ 



(18) 



where {—gf ■ x)~^ represents each non-perturbative gluon bundle, as described in Ref. 
to be exchanged between the quark lines which form the closed loop L[A], its action upon 
L[A] is most conveniently described in terms of a functional cluster decomposition as 



e^A . ^UA] 



exp 
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(19) 



where "connected" requires at least one gluon bundle exchanged between different 
L[^]'s 20|. In this paper we shall be concerned only with the simplest possible applica- 
tion of a single closed loop, and for this we may suppress the linkage operation of (fT9|) . 
while retaining the functional displacement operation of flT5|) . With these simplifications. 



the second line of (fT5l) becomes 

exp {L[A - {gf ■ xV ' (Qi + Q^)]}\a^o = «^P {L[-(^/ " x)'' ■ (Qi + Qn)]}, (20) 
and our eikonal simplifies to 

e«(fe,s) = _^ /"d[^] el/x^ . [det (^/ ■ x)]~^ ■ e^[-(5^->^)"-(2i+2i)l. (21) 

In order to calculate the vacuum loop contribution to (pTI) . we first write a Fradkin 
Representation for L[yl], as in Ref. |3[, 

L[A] = -^Tr /""-e-*'"V(t) /"^H e^/o '^*'-(2/^)-i-^(4)(^(^)) (22) 

^ Jo t J 

X "g-*9/o'^t'<(t')^ni')/d2x^a(xx-xY)A-(x'-^(i')) _ J ^ 

where x|, = (ixg; x^,Xl), h{ti,t2) = 9{ti - ^2)^2 + ^(^2 - ti)ti = ^(ti + t2 - l^i - ^2!), A/'(t) 
is the normalization for the Gaussian functional integral over Va{t'), Tr denotes a trace 
over Dirac and color indices, and the hat notation of a and Q is used to distinguish these 
loop color-variables from those of nucleons I and E. Again, in the interests of simplicity, 
we shall neglect all spin dependence of the quark loop, and, for clarity, have chosen the 
longitudinal and transverse directions of the loop to lie in the respective directions defined 
by the nucleons in their CM. 

With the simplifications of the last two paragraphs, all of the structure that remains 
in our eikonal amplitude arises from that nucleon dependence, Qi and Qi, which has been 
translated into the argument of L in (!2T|) . as its argument A is shifted to — (5'/-x)~^-(Qi+Qn)- 
But this shift occurs in the exponential factor of fl22|) . whose expansion corresponds to 
multiple quark loops exchanged between the nucleons. The simplest, and probably the most 
important effect arises from the exchange of a single quark loop, proportional to the factor Qi 
multiplying Qn, which may be extracted from the quadratic expansion of that exponential 
factor, neglecting tadpole and self-energy corrections to the nucleons. We therefore replace 
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the third hne of ([22]) by 

'+00 „ ft 



+/ j "d.i (^) j "d.2 (^)^''^' / dtit4(tO^^(ii) I dt2t;;,(t2)^'(t2) (23) 
• / d^x'^ a(x_L - x'_l) / d\" a(x_L - x"j_) j d^yf _l a(?/i,± - y'i^^_) j d^t/j _l a(yi,_L - y'^^^) 



■5^'\y'i - sA -x' + v{h)) ■ 6^'\y[ - .2- - x" + v{h)) 
m m 

, db 

Pi/ 



where the generic notation of z'^ = {izq; z'^, zi) for any x'^ or y'^ is understood. Here, x^ 
represents the loop coordinate, which, along with the functional integration of the first 
line of (122|) . will be performed shortly. Notice that various factors of 2, — -i, and g, have 
been combined to produce the coefficient +g'^ multiplying f l23|) . and that this translated 
approximation to L[yl] has become the essence of the desired eikonal function, at least 
before the needed Halpern integrations are performed. 

It will be most convenient to choose the zero of the coordinates yi^ and y-^^ at the 
instant of their CM distance of closest approach, which then corresponds to Si = S2 = 0. 
The argument of each inverse (/ ■ x) is then independent of proper time; consider, e.g., the 
combination {y[ - si^)^ = (yjg - Sif ;?/(_l - Si^,?/i,l - Si^). But in the CM, yi^Q = sif , 
Pj_ ~ 0, and for large momenta ?/l ~ Z/o as pl ~ E, so that {y{ — si^)^ reduces to (0; y[ j_, 0). 
The same argument, with the CM signs of pi^L and j/i^l reversed, holds for the combination 
(z/i-S2^V^(0;yi,^,0). 

This now allows the Si^2 integrals to be performed, and for this a Fourier representation 
of the two delta functions of (l23l) is convenient, which yields 

.Qi'l-{y'l-^')+ikiy'j^-x") _ ^iq-v{ti)+ik-v{t2) 

and where (5(g- Pi) = ^S{qo-qL^) ~ ^S{qo-qi^), S{k-pj) = ^6{ko + ki^S^) ~ ^S{ko + ki^), 
Pl — PiL — ~PnL' ^^^ E = El = E^. Then, the integrations / dxo / dxi^ produce the 
additional factors (27r)^ 6{qo + ko) 6{qi, + /cl), which multiply the previous line, and produce 
a net combination of | ^^-Jj- 6{qo) ^(^0) ^(Ql) (^(^l), so that the remaining q- and k- integrals 
refer to transverse components only. 

Before performing the transverse integrations, it will be convenient to make one further 
simplification, one which appears as a reasonable approximation, but can be justified fol- 

11 



lowing the argument of Appendix B of Ref. ^. This simphfication replaces the arguments 
of each inverse (/ ■ x) factor by their "expected" values yij_ and i/ij,. This step would appear 
to be a reasonable approximation because the a{yi± — y[j_) and a(|/n± — yj±) distributions 
are each peaked about a zero value of their arguments, which essentially forces the primed 
transverse y'j_ coordinates to lie close to their unprimed values. But if appropriate care is 
taken in the evaluation of the x-integrations, one eventually finds the same form of result as 
when this simplification is first performed. Hence, in the interest of clarity and simplicity, 
we now adopt the replacements: x(z/j_) ~^ x{y±)- 

We next evaluate the multiple transverse integrals of (1251) by first writing Fourier trans- 
forms for each of the a{z±) distributions, 

a{z^) = J-^,~a{K)e^-\ (25) 

so that 



(27r)-2 j d'q j d'k j d\'^ j d^x'l j d%^ j d^y;,^ (26) 

■a(x^ - x'^) a(x_L - x" ) a{yi^ - y[^) o(?/i_l - y'^^_) 

.^n-{y'i^-x'J+ik-{y'^^-x'[) 

= (27r)-2|d2g|a(g)na(g)pe'''-« 



where ^{B) will provide a slower fall-off with increasing B than does ip{h) with increasing h. 



Here, B = yix-ynx, and ip 



b) is the modified statement of transverse imprecision introduced 



in [21] and made precise in J], </9(&) — — e~^^^^ , ^ ^ 1. 

The integral of ( l26l) is not the final statement of B dependence, because a term propor- 
tional to Qa qp arising from the evaluation of the functional integral of f l22|) and appearing 
in ( 130|) must still be included. One requires 

M /"d[^;]et/''■(2'^)"'■''^;;^(tl)^^«(t2)5(^^^W)e'^■^'^*'^"'^*'^^ (27) 



which may be accomplished by inserting a Fourier representation of 5'^^\v{t))^ and rewriting 
dZZD as 

Id 5 I d 6 



i dta S gait a) i dtb Sgp{tb) 
■ M [ d[v] e5/^-(2'^)"'-''+^/od*'v(*')[/M{t')+9M{t')] 



(28) 
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where f^(t') = p^ S{t' — t), and 5'^(t') = q^ [S(t' — ti) — S(t' — 12)]. The normahzed, Gaussian 
functional integral of f l28|) is then 

exp |-^ 1^ dt' £ dt" [Ut') + g,{f)] hit', t") [Un + ^,(t")] }, (29) 

and the functional and conventional derivatives of (l28l) . as well as the resulting Gaussian 
J d^p are immediate. Combining all factors, one obtains for the translated and simplified 
L[A] of 022]) the result 

i (^) (1^ /d^g|a(g)na(g)re-^ j^^e~--'n^n^ (3o) 

where we have replaced ti^2 by t ■ 312, and /^ dti_2 by t Jg d^i 2 . 

Eq. (130|) is noteworthy for several reasons, among which is the special way in which the 
manifest gauge invariance of L[y4] is displayed in the automatic cancelation of the quadratic 
divergence associated with the removal of the /^ ^ of (!30|) . In Feynman graph language 
this does not happen automatically, for the divergence of the fermion loop "overpowers" the 
gauge invariance of the basic theory; and one must resort to other measures to remove that 
quadratic divergence. As Schwinger pointed out long ago 8|, in his functional development of 
radiative corrections to QED in terms of proper time variables, such unwanted and improper 
terms never appear in calculations so defined. 

The gauge-invariant divergence of this loop is logarithmic, as expected; and its renormal- 
ization displays the behavior associated with the property of " anti-shielding" , as expected 
in QCD, rather than the "shielding" of QED. This divergence, associated with the lower 
limit of in the t- integral of (1301) . may be described in configuration space by replacing that 
lower limit by a small quantity e, of dimensions of (length)^; in momentum space, this would 
corresponds to a cut-off of A^ = 1/e. It will be convenient to perform the variable change 
t = er, and then rotate contours r — )■ —iz, so that the t- integral of fl30p becomes 

°° d. e-^ [-^-^^^l^-l (^-1^-1)1 n%^ 31) ^\^i2) (31) 



z ' A2 "'' ' A2 



In ( ^ 1 - In 



g2 

l + ;^bl2l(l-bl2l) 
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fi"(0) fi^(O) 



where 312 = 3i — 32, and we have allowed A to become arbitrarily large in the arguments of 
Cl'^ and Cl'^] we have also, for the moment, suppressed the t-dependent integrals coupling d" 
and Cl^ to the Gell-Mann matrices A^. The renormalized coupling of this order g"^ bundle 
diagram may be defined by the relation suggested by flHTl) . as 

A2 



9i{q') = 9' In 



(32) 



gi = g^\n{ — I (33) 



m^ + q^\h2\i^- \h2\). 
which displays the expected QCD form, of an effective, or (partially) renormalized coupling 
that decreases with increasing momentum transfer. And since the g-values expected from 
its subsequent integration are less than the quark mass, and both are understood to be far 
less than any realistic cut-off adopted for A, fl52]) may be most simply approximated by 

A! 

where it is clear that the bare coupling g of the original Lagrangian is smaller than the 
renormalized coupling, in contrast to Abelian QED, where the reverse holds. 

From this example one sees that our formalism is non-perturbative in the sense of sum- 
ming over all gluon exchanges between specified quarks; but that if one of those quark lines 
is part of a closed loop, then a perturbation expansion can be defined involving increasing 
numbers of gluon bundles exchanged between that closed loop and other, specified quarks, 
which may themselves be associated with other quark loops. Can the non-perturbative na- 
ture of our analysis be extended to include all possible L[y4] interactions? We hope to answer 
this very non-trivial question in a subsequent publication. 

The color dependence of fl23l) remains to be treated, and for this it is simplest to return to 
that stage of calculation before renormalization was discussed. There, the factors of Cl'^{ti) 
and Cl'^{t2) remain to be evaluated, which process consists of converting them into Gell-Mann 
matrices A'^ and A'^. It can easily be shown that the commuting factors of Cl'^iti) Cl'^{t2) are to 
be replaced by A'^A'^^(ti— t2)+A'^A'^6'(t2— ti), while, simultaneously, the functional integrations 
over a and Cl have diappeared. 

After renormalization, in which the J dt is effectively evaluated close to its lower limit, 
and where ti^2 =^ t ■ 211,2, as t — > 0, 9(ti — ^2) — > ^(^2 — ^1) — ^ ^(0) = 1/2, and the product 
Cl'^(ti)Cl'^{t2) is replaced by ^{\'^,\'^}. As noted above, for simplicity and ease of presen- 
tation, we have neglected quark spin dependence, and its associated A-dependence, so that 
the Tr operation over both Dirac and color indices yields Tr [^ {A'^, A*^}] = 8 6'^'^. Eq. ( I30ll 
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then reduces to 



^) ^f°* (/• x(yJ)-|::- (/• x(y.J)-'lt (34) 

where the qaQ/s factors of f l5U]) have been replaced by f — ^f- ^§- 
The relevant space-time indices enter here in the form 

(^)(/-x(i)n::-(/-x(n)nt, 05) 

and, remembering the antisymmetry of each element's color and space-time indices, and 
that the a, /3 are transverse indices, Eq. ( 135|) may be rewritten as 

+^{f■ x{i)rXl ■ if ■ xm^'C - ^ (/ ■ ^(i))"Tl ■ (/ ■ xm-tl ' (3^) 

where, because the longitudinal and energy components are far larger than the transverse 
momenta, the /i, u indices correspond to and L only. Using the Minkowski metric, where 
Xla — ^Xoa5 then (/ ■ x)^^\ia = ""■*(/■ x)~^loa- Further, in the CM system, where the lon- 
gitudinal projections of yj and yn point in exactly opposite directions, while the x variables 
depend only upon their respective transverse coordinates; then, the CM longitudinal pro- 
jections of such x^^ will point in opposite directions. In order to have similar, if arbitrary, 
constructions of the x(I) and x(ll), we set X/3l(2) = ~x(2)|^l, so as to bring ( !36l) into the 
form 



(/ ■ x(i))-to ■ (/ ■ ^'(^))"Xl + ^^ ■ ^(^))"'irL ■ (/ ■ xm'Xo \ ■ (37) 



The normalized integrals over J d"x/3L(lI) and J d^x'si^O^), are the same, and are unchanged; 
and since the values of yij_ and |/n± appearing in the arguments of each x serve only to 
indicate that two separate integrations are required, one can interchange those arguments 
in the second term of (1371) to obtain, in place of (1371) . 

- { (/ ■ x{i)rTo ■ (/ ■ xm-'C + (/ • x{i)rXo ■ (/ ■ xW)-'!!} > (^s) 



a result which is explicitly symmetric in a and b and in a and p. 

As in previous discussion of Refs. 2| and [5|], we assume that each x'^ can be represented 



by an angular projection 3"^ multiplying a magnitude R, x'^ = 3^^ R, and we now suppress the 
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result of those normalized angular integrals, assuming that the most significant behavior of 
our results is due to integration over the magnitudes. Of course, such a simplification must 
be checked by detailed, numerical calculation; but this would appear to be a reasonable ap- 
proximation. Note that the index symmetries of f l37p would be enforced by multiplication by 
qaqi3 of (150]) . and by the Q^-Q^, corresponding to color singlet gluon emission and absorption 
of the two nucleons. There is then no difference between the two terms of (l38l) : they are 
both going to give the same contribution, and so fj38|) . after multiplication by the 0"'^, and 
qa, qp is equivalent to 

-2n^Q'q.q, (/■X(I))"X1- (/■^W)~'lo^ (39) 

The attentive reader will notice that there is one aspect of our procedure of obtaining 
an effective potential from an eikonal function which remains to be discussed: what is to be 
done when the eikonal itself contains transverse components of coordinates or corresponding 
momentum transfer? Physically, each component of the initial momentum transfer q^ of 
nucleon I must be transferred to the corresponding component of momentum transfer of 
nucleon E on the other side of the loop, q'i(I) to q'i(ll) and q'2(I) to q'2(ll); in other words, 
a 5ap must appear in fl39|) . either from integrations over the "angular" components of the 
Halpern variables, which we have suppressed, or as a definite statement of our procedure, 
which we now state: All such "free" indices are to be averaged over, a stipulation which has 
consequences in other contexts (renormalization theory and nuclear binding). In the present 
case, it means that qaqp is to be replaced by \q^5ap as is physically necessary. Then, we 
may write the simplified, normalized integrals to be performed as 

N / dRiRl / dR^Rle^^^'^^^^'"^^^ , (40) 

Jq Jo 

where C{B,E) = | (§;] (5^ [(-V|) ■ ^(5)], and, as explained in detail in Ref. U, S^ 

is the scale change needed when passing from the Halpern FI to the individual J d^x'- 

S^ = (fiE)^^. The -Rjj, rather than the -Rjj, result from a factor [Rf^] extracted from 

each determinantal factor [det {gf ■ x)]^ of flT7|) . 

It is the double derivatives with respect to B, the impact parameter between the two 

nucleons, arising from the q'^6a(5 components of the closed-loop integral, which provides the 

sign of a potential that produces nucleon binding; and it is in this qualitative possibility of 

generating a "model deuteron" from two bound nucleons that the possibility of obtaining 

true Nuclear Physics from transversally averaged QCD appears. 
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III. A QUALITATIVE BINDING POTENTIAL 

Before passing to the final steps of the calculation of this potential, we remind the reader 
that this treatment is based on the simplest possible realization of realistic QCD, based on 
a single, massive quark interacting with its complement of SU(3) massless gluons; flavors 
and electroweak interactions, as well as quark and nucleon spins and angular momenta have 
been neglected, for simplicity, and can be added separately, producing definite variations 
of the potential below. The rigorous property of Effective Locality, defined and discussec 
in detail in Refs. |3|, [^ and |5|, immensely simplifies the original Halpern FI of Ref. [2| 
by reducing it, in the present case, to two sets of ordinary integrals; and we have here 
suppressed the "angular" color integrations, retaining dependence only on the magnitudes 
of the reduced Halpern variables, an approximation which must be verified by numerical 
calculation. Nevertheless, it should be of more than passing interest to see just how a 
qualitatively reasonable nucleon potential can appear from such basic QCD. 

Of course, that potential is not meant to suggest that two neutrons will bind, for their 
fermionic nature has been suppressed with the neglect of their spins; nor would it be sug- 
gestive of two protons binding to form a nucleus, because both spin structure and electro- 
dynamics have been omitted. That potential is not yet meant to be compared with precise 
experimental data, except in the sense of its qualitative behavior, producing for two distin- 
guishable nucleons scattering at high relative energies, as well as the possibility of binding 
into a "model deuteron" at lower incident energies. 

With the simplifications and approximations discussed in the preceding Sections, we now 
write flT^ in the form 

/■OO /"OO ^,g c;l 

e^MB^E)^^ di?ii?M dR^Rle"-^''^^''-^-'^^, (41) 

Jo Jo 

where we make the further, simplifying approximation of suppressing the parameter ^ ~ 0.1, 

of (p{b), which was crucial in the construction of quark binding, but would here only slightly 

change the shape of the nucleon binding potential. Setting then (p{b) = (/x^/vr) exp [— (/i6)^], 

one finds 

C{B,E)^^{^y[2^,'B^-].-'^. (42) 

There are several methods of obvious approximation to the integral of (HTi) : 
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1. A change of variables to polar coordinates, Ri = Rsm{6), i?n = Rcos{6), for which the 
radial integral can be done immediately, but the subsequent angular integral requires 
an approximation. 

2. Both i?i.n integrands correspond to a function rising as |-Ri,i| increases from zero, and 
then falling away to as these coordinates become large; and both may be approxi- 
mated by (different) Gaussian approximations. 

3. Both ([1]) and ([2]) lead to rather complicated expressions involving fractional powers 
of complex functions. There is, however, a simpler approximation, available in this 
eikonal context where C{B,E) contains the factor 5^/i^ = {fi/Ey, where E is the CM 
energy of the scattering nucleons. In conventional eikonal representations, there is 
always a dimensionless, energy-dependent, kinematical factor, 7(-E), multiplying an 
impact-parameter-dependent function which is the "true" eikonal function, derived 
from an initial potential function; and if that combination is small, then the final 
eikonal amplitude may have its exponential factor expanded, so that only the linear 
dependence of that exponential is retained. Here, that energy-dependent (5^/i^-factor 
is surely small, but is it the correct 'y{E)7 In Potential Theory and in various forms of 
QFT, the functional form of 'y{E) can vary widely, but we have no precedent here to 
specify the "correct" form of 'y{E) to appear upon the exchange of a pair of gluon bun- 
dles supporting a quark loop. We shall therefore make the simplest choice of adopting 
[lJ,S{E)]^ as our tentative 7(-E); and at the very end of the calculation return to see if 
this choice is consistent with the order-of-magnitude of our qualitative potential. 



We now expand to first order both the eikonal amplitude of flT3|) . which is the left-hand- 
side of (HTj) . and the exponential factor containing C{B, E) of (HT1) . so that 



iX{B,E) = iX^{E) = -iC{B,E)Afj\ (43) 



J= I di?i?3e3^, 
'o 

•4' 



AfJ^ = -I- 



where the "true" eikonal function is 



M /^R 



^"(e Ut)['-^'^']^"^- (^^) 
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This situation differs from that of the quark binding calculation of Ref. [5|, where the large 
impact parameter of interest generated a large eikonal function, but a small amplitude; here, 
both the eikonal and the amplitude are small. 

The relation between the eikonal function and the effective potential is 



X=-/ dz^V{B + P^zj^) (45) 

and as explained in Ref. p], the eikonal is real for a purely scattering potential, V = Vs; 
but for a potential which can lead to binding, or to the production of other particles, the 
potential chosen must have the form V^ = Vs — iVe, so that the eikonal which corresponds 
to binding is imaginary. The reason is unitarity, since if extra particles, or a new bound 
state can be produced, the amplitude of the initial state must be reduced. In our reversed 
situation, starting from the construction of a QCD amplitude, we find a clear signal of a 
binding potential, with Vb appearing as a real quantity, 

/+00 1 2 

dzL ^B = Y^ If (2 - f^'B') e~^'^''\ (46) 

To obtain the effective potential one first calculates the two-dimensional Fourier transform 
of —iX{B), which can be expressed as proportional to 

flpl kl U^Be^^-^~>'"^"'\ (47) 

then continue k^ to three dimensions, k^ — )■ k\ + k^^, and calculate the three-dimensional 
Fourier transform of f H7|) . which yields, after removing the factor 7(-E) = (/i(5)^, 

V{r) ~ c (^ /i [2 - A'] e-"^, (48) 

with c = I (27r)^'^/^. At high energies and large momentum transfers, this potential when 
multiplied by (— i) corresponds to an effective scattering potential. 

The form of this potential is sketched in Fig. [3l and it will look familiar to those who have 
inferred a nucleon potential from experimental scattering data, starting with the potentials 
of the 1951 paper of Jastrow [12]. It must be noted that this potential is not meant to 
be relevant at distances /ir < 1, which is where the multiple gluon exchanges of the gluon 
bundles of Fig. [21 as well as those of the omitted gluon-binding interactions of each triad of 
bound nucleons take place. And of course, we have neglected electromagnetic effects, as well 
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FIG. 1: A gluon bundle exchanged between two nucleons 




FIG. 2: Quark loop exchange through gluon bundles between two nucleons 



as all spin and angular momentum modifications, which can be included in more detailed 
estimations. 

We have two parameters at our disposal, the mass scale /i ~ m.,^, and g^/Air, which can 
be chosen so as to produce a ground state with a binding energy of -2.2 MeV. Of course, 
from the crudeness of the approximations made in our various estimations, we would be 
happy to obtain a binding energy to within a factor of 10 of this numerical result, but as it 
happens, we shall do somewhat better. The corresponding calculation is demonstrated in the 
next Section, using the elementary Quantic technique [9(] of estimating a ground state. But, 
simplifying approximations aside, this is clearly a potential which can bind a pair of distinct, 
uncharged nucleons; and it is obtained analytically, from basic transversally-averaged QCD. 



FIG. 3: Nucleon potential 
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IV. BINDING ESTIMATIONS 

The non-relativistic Hamiltonian of two, equal mass particles interacting with the above 
potential is 

E = ^ + V{r)^^ + V{r). (49) 



One can write this non-relativistic energy in dimensionless form as 

^ + ^[2-y2]e-^, y = ;.r, (50) 



E{y) _ 1 Vo ro ..21 
M y^ M 



where we have set /i = m.^. In units of MeV, M = fi'^/m = 18.2, g^/4:7T = K ■ 10, and, 
combining all the relevant factors of the previous paragraphs, Vq = 14.3 K (MeV). 

The Quantic method |_9'] of estimating a ground state is to rewrite p as 1/r = fx/y, to 
find the minimum of E{y), and use that minimum point yo to define E{yo), which is to be 
interpreted as a qualitative estimate of the ground-state energy. The minimization statement 
is given by the vanishing of the derivative of (l50l) at y = yo, 



= -^ + T7l/obo-4]e-^, (51) 



yl M^' 

and the customary way of solving such a problem is to solve (ISTl) for yo, and then substitute 
that value of yo back into ( l50l) to obtain the binding energy. But since K is essentially 
unknown - one might guess it to be on the order of 1, representing a strong, nuclear force 
- and because we do want to represent the bound state energy as E{yo) = —2.2 Mev, let us 
use that number together with the value of M to solve for yo; and then solve for the value 
of K. 

To do this, combine ( l50l) and ( 15T1) in such a manner that the exp [—yo] factors of both 
equations are canceled, which produces 

1^1 



Vo (2/0 - 4) 



^ + ^0M 



2(^0-4) (52) 



which is a cubic equation in y"^. From the graph of Fig. [3l one sees that the minimum of 
the potential lies close to yo = 2, which suggests that the minimum of the energy should 
be somewhat larger; and this suggests the choice yo = 3 + A as a reasonable choice for 
the approximate solution of ( l52l) . retaining terms of no higher order than A (under the 
subsequently verified assumption that A^ ^ 6|A|). This leads to the result: A ~ .33 and 
yo ~ 3.33. Upon substituting this value of yo into (|50|) there follows K ~ 1.25, which provides 
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the expected order of magnitude for a strong-coupling process. One may expect that when 
contributions from quark and nucleon spins are included, that number will decrease slightly, 
retaining its strong-coupling character. 

V. SUMMARY AND SPECULATION 

While the arguments put forth above are concerned with a realistic version of QCD, 
and have for simplicity neglected flavors, and electroweak interactions, quark and nucleon 
spin dependence, and have suppressed several "angular" integrations, all of which can be 
restored, as desired, the result is an explicit, model "deuteron" potential, of sufficiently short 
range and of the right order of magnitude to be considered as a qualitative derivation of 
nucleon-nucleon forces from basic realistic QCD. 

Our tentative choice oi'y{E) = (/i 5)^ = (fi/E)^ has turned out to be qualitatively correct; 
and, in an eikonal context, this is interesting because it suggests that for the exchange of a 
composite object - in this case, the gluon bundles supporting a quark loop - between two 
"scalar" nucleons, the 7(-E) factor is not just what one typically finds when exchanging scalar 
quanta, {m/Ey, but retains the memory and has a signature of the composite structure 
being exchanged, 'y{E) = {fi/m)'^ {m/ E^ . 

The above analysis should be almost immediately applicable to high-energy nucleon- 
nucleon scattering; and it will be interesting to see if the result of that calculation corresponds 



to the physical arguments recently suggested by Islam [ll|. 

Generalizations of this two-nucleon deuteron model to the construction of heavier, stable 
nuclei may well be possible, and might provide at least a partial basis for the nuclear shell 
model and the independent boson (IBM) model. In the first case, one would ask how many 
nucleon-generated gluon bundles can be attached to a single quark loop; and for the IBM 
model, asking how effective would attractive pairwise interactions of the deuteron form be 
when exchanged between nucleons in a three-dimensional array. 

Finally, on a more fundamental level, it will be most interesting to see just how the struc- 
ture of renormalization theory turns out for realistic QCD, the theory which has, built-in, 
quark transverse imprecision. From the experience gained in our work so far, the simplifi- 
cations in which non-perturbative gluon exchanges organize themselves into gluon bundle 
exchange displaying Effective Locality suggest that truly non-perturbative renormalization 
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will turn out to be simpler than that of QED. We hope to answer this question in the near 
future. 
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